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Abstract 

Based on the Vlasov-Maxwell equations describing the self-consistent nonhnear 
beam dynamics and collective processes, the evolution of an intense sheet beam 
propagating through a periodic focusing field has been studied. It has been shown 
that in the case of a beam with uniform phase space density the Vlasov-Maxwell 
equations can be replaced exactly by the hydrodynamic equations with a triple 
adiabatic pressure law coupled to the Maxwell equations. We further demon- 
strate that starting from the system of hydrodynamic and Maxwell equations 
a set of coupled nonlinear Schrodinger equations for the slowly varying ampli- 
tudes of density waves can be derived. In the case where a parametric resonance 
between a certain mode of density waves and the external focusing occurs, the 
slow evolution of the resonant amplitudes in the cold-beam limit is shown to 
be governed by a system of coupled Gross-Pitaevskii equations. Properties of 
the nonlinear Schrodinger equation as well as properties of the Gross-Pitaevskii 
equation are discussed, together with soliton and condensate formation in intense 
particle beams. 



1 Introduction 



One of the main goals in the commissioning and operation of modern high energy acceler- 
ators and storage rings is the achievement of higher and higher beam currents and charge 
densities. That is why, of particular importance are the effects of intense self-fields due to 
space charge and current, influencing the beam propagation, its stability and transport prop- 
erties. In general, a complete description of collective processes in intense charged particle 
beams is provided by the Vlasov- Maxwell equations for the self-consistent evolution of the 
beam distribution function and the electromagnetic fields. Although the analytical basis 
(as mentioned above) for modelling the dynamics and behaviour of space-charge dominated 
beams is well established, a thorough and satisfactory understanding of collective processes, 
detailed equilibrium and formation of patterns and coherent structures is far from being 
complete. 

The present paper may be regarded as the first (to our knowledge) attempt to take 
a view at the description of the evolution and the collective behaviour of intense charged 
particle beams from an entirely different perspective, as compared to the ones available in 
the literature. We will be mainly interested in describing the slow evolution of some coarse- 
grained quantities that are easily measurable, such as the amplitudes of density waves. Due 
to the nonlinear wave interaction contingent on the nonlinear coupling between the Vlasov 
and Maxwell equations, one can expect a formation of nontrivial coherent structure that 
might be fairly stable in space and time. Here, we show that solitary wave patterns in 
the beam density distribution are an irrevocable feature, characteristic of intense beams. 
Moreover, density condensates in the special case where a parametric resonance in terms 
of wave frequency between a particular mode of the fundamental density waves and the 
external focusing occurs, can be formed. 

The paper is organized as follows. It was previously shown pQ that in the case of a sheet 
beam with constant phase-space density the Vlasov-Maxwell equations are fully equivalent 
to a hydrodynamic model with zero heat flow and triple-adiabatic equation-of-state. For 
consistency, in section 2, we repeat the derivation from a slightly different perspective and 
then use the hydrodynamic equations as a fundament for the subsequent analysis. In section 
3, we consider the case of a smooth focusing where the time variation of the /5-function (re- 
spectively, the time variation of the density envelope function) can be neglected. Further, we 
demonstrate that starting from the hydrodynamic equations, and using the renormalization 
group (RG) technique [21 El IH E] a system of coupled nonlinear Schrodinger equations for 
the slowly varying amplitudes of density waves can be derived. The purpose of section 4 is 
twofold. First, we study the case where a parametric resonance between a particular mode 
of density waves and a resonant Fourier harmonic of the external focusing (a resonant har- 
monic from the -y/^-function Fourier decomposition) occurs. It is shown that under certain 
conditions the evolution of the resonant amplitudes of the forward and backward density 
waves is governed by a system of coupled Gross-Pitaevskii equations (see e.g. the review ^ 
and the references therein). Secondly, it is demonstrated that in the non- resonant case the 
results obtained in the smooth focusing approximation coincide with the ones for a periodic 
focusing up to second order in the formal perturbation parameter. In section 5, the reduction 
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of the infinite system of coupled nonlinear Schrodinger equations to a system of two coupled 
nonlinear Schrodinger equations is discussed. Such reduction is possible if one neglects the 
effect of all other modes and takes into account only the self-interaction of a single mode 
with a particular wave number. In the non-resonant case it is shown that the renormalized 
solution for the beam density describes the process of formation of holes (cavitons) in intense 
particle beams. The case where a parametric resonance is present is more interesting. In 
the cold-beam limit it is demonstrated that the evolution of the forward and the backward 
resonant wave amplitudes can be well described by a system of two coupled Gross-Pit aevskii 
equations. Finally, conclusions are drawn in section 6. 



2 Derivation of the Hydrodynamic Model 

We begin with the Hamiltonian describing the one-dimensional betatron motion in the pres- 
ence of space-charge field 

H = ^p^ + + V{x; 9) + ^fsciS; 9), (2.1) 

where e is the electron charge, R is the mean radius of the accelerator, G{9) is the focusing 
strength of the linear machine lattice, V{x; 9) is a contribution coming from nonlinear lattice 
elements (sextupoles, octupoles, etc.) and Eg and Ps are the energy and the relative velocity 
of the synchronous particle respectively. In addition, ipsc{x; 9) is the self-field potential due 
to space-charge, which satisfies the one-dimensional Poisson equation 



d'^^Psc _ _en 



dpf{x,p-9). (2.2) 



Here n = Np/Vt is the density of beam particles [Np is the number of particles in the beam, 
while Vt is the volume of the area occupied by the beam in the transverse direction), eo is the 
dielectric susceptibility and f{x,p]9) is the distribution function in phase space. Equations 
()2.H) and ()2.2|) can be written in an alternate form as 

H = ^f + + V{x- 9) + A^(2; 0), (2.3) 



dx"^ 



- J dpf{x,p;9), (2.4) 



where 



A = -1^, (2.5) 



Re^n 

is the beam perveance, and ipsc = enip/eo. Note that the beam perveance A is a dimensionless 
quantity. 

Next, we perform a canonical transformation 

x = x^, p=^{p-ax), (2.6) 
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defined by the generating function 

where a and (5 (and 7) are the well-known Twiss parameters satisfying the equations 

As a result, we obtain the new Hamiltonian 

^ = + x^) + V{x- 9) + \^U{x- 9), (2.9) 

where the self-field potential U (x; 9) ((p = y/^U^ satisfies the equation 

d'^U r 

-g^ = -J dpf{x,p-9), (2.10) 

and 

is the derivative of the phase advance with respect to 9. 

We are now ready to write the Vlasov equation for the one-particle distribution function 
f{x,p; 9) in the two-dimensional phase space (a;,p). It reads as 

^ + _ f + ^ + XS^A ^ = 0, (2.12) 

89 dx \ dx ^ dx ) dp ' 

and should be solved self-consistently with the Poisson equation ()2.10|) . Following Davidson 
et al. PP, we consider the case where the distribution function f{x,p]9) is constant (in- 
dependent of X, p and 9) inside a region in phase space confined by the simply connected 
boundary curves p(+)(x; 9) and P(-)(x; 9), and zero outside. In other words, 

f{x,p-9)=C, for p(_)(x;^) <p<p(+)(x;^), (2.13) 

and f{x,p;9) = otherwise. The proof that such solution to the Vlasov equation (j2.12p 
exists can be found in appendix A. 

The equations for the boundary curves P(±){x] 9) can be derived as follows. We substitute 
the explicit form 

/(x,p; 9)=c[n{p- pi-){x- 9)) - n{p - p(+)(x; 9))] (2.14) 

of the distribution function into the Vlasov equation (|2.12|) . Here ?i(-z) is the well-known 
Heaviside function. Using the fact that the derivative of the Heaviside function with respect 
to its argument equals the Dirac (5-function, we obtain 



dx 



4 



dV 



\ ox ^ ox 



0. 



[2.15) 



Multiplying equation ()2.15p first by 1 and then by p, and integrating the result over p, we 
readily obtain 

+ f -pU) = 0, 



ld_ 

2de' 



de 

{ph - pU) + l^fe - pU) = -(m) - Pi-)) + ^ + ^V^^ ) ' (2-17) 



(2.16) 



dV 



dx 



dx I 



8nj_ 

dx"^ 



P{+)-P{-)h 



:2.i8i 



It is convenient to recast equations (j2.16j) - (j2.18|) in a more familiar form, widely used 
in hydrodynamics. Let us define 



(2.19) 



C 



pv = dppf{x,p;9) = -{pl^^-pf_^) 



(2.20) 



where p{x] 9) and v{x;6) are the density and the current velocity, respectively. From the two 
definitions above, it follows that 



v = l{p{+)+P{-))- 



[2.21) 



Defining further the pressure V{x;6) and the heat fiow Q{x;9) and using equation ()2.2ip . 
we have 



V 



(lp{p-vff{x,p] 9) 



C{p(+)-P(~ 



12C2 

oo 

Q= j dp{p~v)^f{x,p;9) = 0. 



(2.221 



(2.23) 



From equation (|2.19|) . it follows that the pressure can be expressed as 

1 



Po 

In addition, it is straightforward to verify that 



PI 12C2- 



(2.24) 



(2.25) 



which provides a closure for the equations (j2.16|) and (|2.17|) governing the evolution of the 
moments (boundary curves). With all the above definitions and relations in hand, we can 
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write the system of hydro dynamic equations (supplemented by the Poisson equation) in the 
form 

|h-4(H = 0. (2.26) 
|,H + 4(^> + ,.^)^-p(.. + f H-Av^g, (2.27) 



-3^ - 

Using the continuity equation ()2.26|) . we cast the system of equations (j2.26|) - ()2.28p in its 
final form 



| + |(P.)^0. (2.29) 

ft. ^ ^ dv ^ 3 , 2\ ^ 13 av ^ti^dU 

dr dx ^ dx^ ' R dx R dx' 

-P, (2.31) 



dx'^ 

which will be the starting point for the subsequent analysis. Here 

2pg 



(2.32) 



is the normalized thermal speed-squared, and r = x(0) + tq is the new independent "time"- 
variable. 



3 Renormalization Group Reduction of the Hydrody- 
namic Equations 

Before proceeding further, we make an important remark. Let F{d) be a periodic function 
of 6 with period 27r. Noting that the phase advance x(6') can be represented in the form 

x{e) = ye + Xp{e)~xp{T^). (s.i) 

where v is the betatron tune, and Xp(^ + 27r) = Xp(^)) we can expand F{6) regarded as a 
function of x (respectively r) in a Fourier series in the new variable x (respectively r) as 
follows 

F{e)= f; Xexpl'^n^y (3.2) 
The expansion coefficients An are given by the expressions 

An = ^ I dxF(0) exp f-m^y (3.3) 
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Using the definition of the phase advance (|3.1|) . and choosing tq = Xp(7r), we can rewrite the 
Fourier expansion (|3.2|) as 

r 



F{e)= ^nexp 



m- 



(3.4) 



where the expansion coefficients i3„ are given by the expressions 



R 

27ny 



ciU ^^^^ e exp 



(3.5) 



In what follows, we consider the case where the external potential V is zero, so that the 
set of equations to be analyzed acquires the form 



9p d , , ^ 

dv dv 2 ^ ( 2\ ■> / 

f)-T f)-r \ ) ^ Qx 



dr dx dx 



where 



dx"^ 

CO , ^ 

9{t) = «n.exp (in- 

n=— oo ^ 

1 ^ 

in = 7^ f d^7^e-^"^xp 



-in 



(3.6) 

(3.7) 
(3.8) 

(3.9) 

(3.10) 



Further, we assume that there exists a nontrivial solution to equations (j3.6|) - ()3.8|) in 
the interval x G a;(+)^ , and that the sheet beam density is zero {g = 0) outside of 

the interval. We introduce also the ansatz 

1 X x^ 

p{x;t) = - + eR{x;T), v{x;t) = - —+eu{x;r), U{x;r) = - —+eU{x]T), (3.11) 

where the envelope function S{t) is a solution to the equation 



dr2 



+ £ = \giT), 



and therefore can be represented in explicit form as 
£(r) = Aao + A^6„e*"^/^ 



1 — V? jv"^ 



(3.12) 



(3.13) 



It enables us to rewrite equations ()3.6p - ()3.8p as follows 



+ eu— + ev^T—(R^)=-\g—, 



dr ~^ £ dx ~^ £ dx 



(3.15) 



-R. 



Next, we differentiate the first equation ()3.14|) of the above system with respect to r and the 
second equation ()3.15p with respect to x. After summing these up, we obtain 



- 2v. 



T 



d^R 
dx"^ 



+ Xg£R 



£ 



d 



+ e£ 



dx^ \ 2 



+ vlR'^ - e£ 



dx"^ ^ ^ dr 

Equation ()3.16|) supplemented with equation ()3.14p rewritten as 

^dR du ■ d , ^ 9 , 

^1^ + ^irixR) + t£—{Ru) = 0, 
OT OX ox ox 



d_ 



'£liRu) 



(3.16) 



(3.17) 



comprises the basic system of equations for the analysis in the subsequent exposition. 

In what follows, we consider the case of smooth focusing, where the time variation of 
g{T) [£{t), respectively] can be neglected. In the next section, we will show that when the 
frequencies of all fundamental modes are sufficiently far from a parametric resonance, this 
assumption holds true to second order in the perturbation parameter e even in the case when 
such time variation is present. Thus, the basic equations ()3.16|1 and ()3.17|1 acquire the form 



2vi 



d^R 

dx"^ 



+ £qR — €£q 



dx"^ 



u 



+ vlR^ - e£n 



drdx 



(Ru) 



where 



^ dR du ^ d , ^ ^ 
or Ox Ox 



£o — XaQ. 



(3.18) 

(3.19) 
(3.20) 



Let us further assume that the actual dependence of R and u on the independent variables 
is given by 

R = Rix,^;T), n = n(x,e;r), (3.21) 

where 

^ = ex (3.22) 

is a slow spatial variable. Therefore, the basic equations (|3.18|) and (j3.19p can be rewritten 

as 



.2^ 



^""hdx^ 



+ 2e 



8 



u 



2 d2 



^ OR du du ^ ( d d \ ,^ 



° yOrdx ' drd^ 
0, 



\{Ru) 



(3.23) 



(3.24) 



Following the basic idea of the RG method, we represent the solution to equations ()3.23|) 
and ()3.24|1 in the form of a standard perturbation expansion [Zj in the formal small parameter 

e as 

oo oo 

R=Y.e'^Rk, u = J2e''u,, (3.25) 

k=0 k=0 

The zeroth-order equations ()3.23|1 and ()3.24|1 read as 



\2„2 



Rq „ 9 R{ 



2vi 



° + X^alRo = 0, 



dRo duo „ 

Aao^:^ h ^r— = 0. 

ar ox 

Their solutions can be found in a straightforward manner to be 



Uo{x,i;T) 



2tT , UJm 

T2 2^ T;r 



a 



m 



.(-) 



(3.26) 
(3.27) 

(3.28) 
(3.29) 



Here and are constant complex amplitudes (to this end dependent on the slow variable 
^ only) of the backward and the forward wave solution to equation (j3.26|) . respectively. These 
will be the subject of renormalization at the final step of the renormalization group procedure 
resulting in RG equations governing their slow evolution. Furthermore, 



z^^'> (x; r) = ± max, 



27r 



Aao 



27r 



The discrete mode frequencies Um are determined from the dispersion relation 



2 _4 



1 + JsTW, 



K 



'2^ 



It can be easily verified that the above choice of parameters leads to 



'-(+) 



dxRnix; 



(3.30) 



(3.31) 



(3.32) 



which means that linear perturbation to the uniform density Pq = £q^ average to zero and do 
not affect the normalization properties on the interval {x'^~\x^^^^. In addition, the following 
conventions and notations 



-UJr. 



A. 



A* 



B_ 



B* 



(3.33) 
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have been introduced. 

The first order perturbation equation for Ri reads as 



+ Ri 



- E 

<7 -An 



E 



7, 



where 

The solution for Ri is readily obtained to be 



uJm + ^n){ — ± — ] + {m±n){K'± 



-L R B e 

^ /mn-'-'m-L'ne 



^2 I '^Tn^n 



(3.34) 
(3.35) 



i^i(x,e;r) = ^E-f%^e-^"-§^e-^'' 



- E 



where 



(±) _ /mn 



Note that the (infinite dimensional) matrices a^^^ are symmetric, i.e. 



mn nm ' 



a^^^ - 



(3.36) 

(3.37) 
(3.38) 

(3.39) 



Having determined Ri {x; r) , the first-order current velocity Ui {x; r) can be found in a 
straightforward manner. The result is 



Ui{x,C;t) 



E 



2m 



1 / dAm. ,-,(+) dBm^ ■{-) 



^0 ^^^^ \ 



27r2 „ 



where 



/3(±) _ _1_ I (±) 

m n m±n 



= 



(3.40) 
(3.41) 
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Note that the (infinite dimensional) matrix possesses the same symmetry properties as 
those displayed by equation (|3.39p possessed by the matrix a^~^\ while is antisymmetric 
(and evidently = 0). 

In second order, the equation for R2{x,^;t) acquires a form similar to that of equation 
fl3.34p . It is important to emphasize that this assertion holds true in every subsequent or- 
der. Each entry on the right-hand-sides of the corresponding equations can be calculated 
explicitly utilizing the already determined quantities from the previous orders. The right- 
hand-side of the equation for R2{x,^; r) contains terms which yield oscillating contributions 
with constant amplitudes to the solution for R2{x,C,',t). Apart from these, there are reso- 
nant terms (proportional to e*^™ (^;^)) leading to a secular contribution. To complete the 
renormalization group reduction of the hydrodynamic equations, we select these particular 
resonant second-order terms on the right-hand-side of the equation determining R2{x,^] r). 
The latter can be written as 



+ 



2iK^T ^ , (+) . d^B^, , J-) 



m iz^ 

e " 



,{-) 



where 



r(+) 



m 



r(-) 



m 



m n 

UJm. UJr, 



mn 



A 



(-) 

mn 



m n 



mn 



Some straightforward algebra yields the solution for i?2(x,^; r) in the form 



^2{x,i\T) 



X ^ III 

h< \ de 



+ 



(3.42) 

(3.43) 
(3.44) 



1 id'' Am . 



+ 



,(-) 



(3.45) 

where non-secular oscillating terms have not been written out in full explicitly. The final 
step is to collect the terms proportional to the fundamental modes e*^™ and e*^™ in all 
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orders and renormalize the amplitudes and B^. As a result one obtains the following 
RG equations 



2iuj. 



OA. 

dr 

dBrr 



— 2im- 



+ —^Am'^ (Tml\A\ +rm2|-Bn|), (3.46) 



dr 



(3.47) 



4 The Parametric Resonance 



Let us now address the system of equations (j3.16p and (j3.17p . Without loss of generality, we 
assume that the time variation of g{T) can be treated as a second-order perturbation (which 
is usually the case), that is 

g{T) = ao + e 2^ a„e ' . 
The same holds true for the envelope function S{t) 



(4.1) 



£(r) = Aao + e'A5:6„e^^^/^ 



(4.2) 



Thus, the basic equations ()3.16|) and ()3.17p can be rewritten in the form 



d 



dx"^ 



\xu) 



dT\ dr ^ 



-2vi 



dx"^ 



+ Xg£R 



dr 

OR du 



OT OX OX OX 



0, 



(4.3) 



(4.4) 



Let us reiterate that the assumption concerning the second order of magnitude in e of the 
time variation of g{T) and £{t) does not restrict the generality of the subsequent analysis. 
If this variation were of first order in e, the proper perturbation parameter to use would be 
s/e instead of e. In addition, the variables R and u have to be rescaled accordingly 



R 



R_ 

7^' 



u 



u 

7^' 



(4.5) 



In terms of the ansatz 1)3.211) and ()3.22|1 . equations ()4.3j) and ()4.4|1 become 



- 2u 



T 



dx' 



+ 2e 



92 



92 



\R + \g8R 



dx^ 



2t 



92 



d_ 



£2 



d 



dx 



12 



u 



+ 1 - e-^ 



_9_ 



du 



ar ax 



d 



9. + ^af]<^''> 



0. 



(4.6) 



(4.7) 



Although the general case can be in principle treated through more labour-intensive manip- 
ulations, for the sake of simplicity in what follows, we select a particular mode with mode 
number m. The zeroth-order solution can be written in the form 



R,{x,i-T) = A(Oe^^'^'(^'-) + 5(0e^^^"'(^'^) + C.C., 



2ttuj 



C.C., 



(4.8) 
(4.9) 



where again 



z^'^\x, t) — UJT ± max, 



= 1 + K'^w?. 



The particular mode is chosen such that an exact parametric resonance (if possible) 

2uj^— (4.10) 



no 



occurs for some integer no. The first order perturbation equation for Ri can be written in 
the form 



9t2 



dx^ 



+ Ri^ 



2imK^ (dA ■ 



:(+) 



- (4a;^ - l) (AV-'"' + B'e'-'-') + ^AB*e<^'^'-^'-') + c.c. (4.11) 



^(4a;^-l)(AV-'^'+SV-^-^) + ^ 

Here Y^' implies exclusion of the harmonic with n = Uq from the sum. It is important to 
note that terms giving rise to secular contribution due to the parametric resonance vanish 
identically in the first order. The solution for Ri can be written explicitly as 



Straightforward calculations yield the solution for u\ 
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iXrir 



V 



.(-) 



iuUiu' - h/2){A'e''''^' - B'e'''' ' ) + c.c 



(4.13) 

In second order, we take into account terms providing secular contribution to the solution 
only. Thus, we write the equation for R2 as 



+ R2 = —t\ -^e'^ + -^e'' 



TC 



+ 



a2 I 9^ 



9e 



2 2 



.(-) 



3^ 



+ 



167r^ 



\B\^Ae'''*^ + |A|'Se'"*"') +C.C.. 



(t2(4cu2-3) 

Similar to equation (j3.45p . we obtain in a straightforward manner 



(4.14) 



K^m^ Jd^A d^B 



+ 



2zcj3a2 I 5^2 



9f 



;4-i).v-.(4-i).v- 



(-) 



(l6a;^ - lloo^ + l){\AfAe'''''' + \BfBe"''' 



2iuj 3o"2 

T 167r2 



2zcj(t2(4^2 _ 3) 



|5|2^e-(+' + \AfBe''^'^)+c.c. 



(4.15) 
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Collecting once again the zeroth-order term together with all secular terms in higher orders 
and renormalizing the amplitudes A and B, we obtain the RG equations 



^. dA K^dA d^A Xa 

2ujj— 2im — 



dr a dx a'^uP' dx^ 



^^'^'-m2c.2|6„jVA - —A\^u' - IW + l)\AfA + ^/^f \B\^A, (4.16) 



dB K^dB d^B Ad, / .\( d \ 

OT a ox a^io^ ax^ vr ^ ^ \ ax j 

+ ^mV|6„J VS - ^^(l6u' - IW + l) \BfB + ,,ff ^. (4.17) 

5 The Nonlinear Schrodinger Equation for a Single 
Mode 

Equations ()3.4(i|l and ()3.47j) represent a system of coupled nonlinear Schrodinger equations 
for the mode amplitudes Am and Bm- Neglecting the contribution from modes with n ^ ±m 
and introducing a new amplitude Bm{x] r) = Bm{—x] r), for the single mode amplitudes Am 
and Bm, we obtain the equations 

.dB d'^B Sn^ 



dr dC cj(t2(4cj2 - 3) 



{\Af -G\Bf)B = 0, (5.2) 



where 



C = V2uj{ mKr + , G = ^^^^^ (l6cu^ - 1 Icu^ + l) , (5.3) 



the index m and the tilde sign over the new amplitude B has been omitted. Clearly enough, 
equations (|5.1|) and (j5.2|) follow directly from the system (j4.16|) - (j4.17p if 6„q is set to 
zero. This implies that in the case where a parametric resonance does not occur in the 
original system, the smooth approximation is valid to second order in the formal perturbation 
parameter. Moreover, it is straightforward to verify that G is always positive. The system 
of two coupled nonlinear Schrodinger equations ()5.ip and ()5.2|) is in general non integrable. 
Its integrability is proven by Manakov |HI only in the simplest case of G = — 1. 

If one of the amplitudes {B or A) in its capacity of being a particular solution to the 
corresponding nonlinear Schrodinger equation is identically zero, the equation for the other 
amplitude (say A) becomes 

- |1 + n^l^A = 0. (5,4) 



15 



where 



47r^ 



(5.5) 



In nonlinear optics equation ()5.4|1 is known to describe the formation and evolution of the 
so-called dark solitons [Hj. In the case of charged particle beams these correspond to the 
formation of holes or cavitons in the beam. The solution to equation ()5.4|) can be written as 



r(0 

7f' 



,i[nr+e(0] 



where 



for all c and 



n — 2a^sech^(a^), 



0{^) = arctan 



cr, 



2a , , . ■ 
— tanh(a^) 



1 



V2n- 



(5.6) 

(5.7) 
(5.8) 



provided n > 

Let us now examine equations ()4.16|) and ()4.17|) . In the cold-beam limit f t — > (or 
equivalently u ^ 1) the second term on their right-hand-sides can be neglected as compared 
to the other terms. Therefore, we can write 



dA 



+ 



■^^ 2 217, |2 2 7 

— m uj p„ , X A 



A 



A + 



a2(4cu2-3) 



B 



A. 



(5.9) 
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Or (j'^ijp- dx^ 



(l6cj^ - llcu^ + 1 



B 



B + 



^"0 I 

167r2 



a2(4cj2-3) 



5, 



where the local gauge transformation 



A = A exp 



imuj 



-T + aujx 



B = Bexp 



fmK^ ^ 

imoj r — auJx 



(5.10) 



(5.11) 



has been performed on the amplitudes A and B. Equations ()5.9|) and ()5.10|) represent a sys- 
tem of two coupled Gross-Pitaevskii equations, known to govern the formation of condensate 
structures in atomic gases confined in magnetic traps. In addition, superfluidity in helium 
as well as patterns in the gas of paraexcitons in semiconductors are considered as a possible 
manifestation of Bose-Einstein condensation. It is remarkable that under certain conditions 
similar phenomenon can be observed in space-charge dominated beams. 
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6 Concluding Remarks 



The analysis performed in the present paper is based on the Vlasov- Maxwell equations for 
the self-consistent evolution of the beam distribution function and the electromagnetic fields. 
We considered the propagation of an intense beam through a periodic focusing lattice. It 
has been proven that a special class of solutions to the Vlasov-Maxwell equations exists, 
which defines a uniform phase-space density inside of simply connected boundary curves. 
These solutions provide an exact closure to the hierarchy of moments resulting in a set 
of hydrodynamic equations for the beam density (continuity equation) and for the current 
velocity. The latter is characterized by a triple adiabatic pressure law. The major drawback 
of the model discussed here is that due to the constancy of the distribution function inside 
of the evolving region of phase space, it fails to take into account an important feature such 
as the well-known Landau damping. Inclusion of the Landau damping in the description can 
be achieved correctly only by direct analysis of the Vlasov-Maxwell system [TU| . 

Further, we studied first the case when the smooth focusing approximation applies. Based 
on the RG method, a system of coupled nonlinear Schrodinger equations has been derived for 
the slowly varying amplitudes of interacting beam-density waves. Under the approximation 
of an isolated mode neglecting the effect of the rest of the modes, this system reduces to 
two coupled nonlinear Schrodinger equations for the amplitudes of the forward and of the 
backward wave. The particular solution to the latter system asserting that the amplitude 
of either of the waves (the forward or the backward) can vanish identically, leads to a single 
nonlinear Schrodinger equation for the other amplitude. The latter is characterized by a 
repulsive nonlinearity and therefore, it describes the evolution of holes in intense charged 
particle beams. 

The analysis of periodic focusing clearly showed that the results obtained in the case of 
smooth focusing remain unchanged up to second order in the perturbation parameter if a 
parametric resonance between a particular mode and an appropriate Fourier harmonic of the 
external focusing does not occur. If however, an exact parametric resonance takes place, it 
was shown that the evolution of the wave amplitudes of the resonant mode in the cold-beam 
limit is described by a system of coupled Gross-Pitaevskii equations. Quite remarkably, it 
was found that there exist a possibility of formation of density condensates in space-charge 
dominated beams. 
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Appendix 

In order to prove that the expression ()2.14|) represents the solution of the Vlasov equation 
()2.12|) . we write the equations for the characteristics of the latter in the form 

1 Xp \ ox ^ ox J U 

Let us also assume that the distribution function f{x,p; 6q) at some initial time 60 is given 
by 

f{x,p;eo) = fo{x,p)=C[n{p~pf}^{x))-n{p-p^l^^{x))]. (A.2) 

Suppose now that we have been able to solve the equations for the characteristics ()A.1|) . and 
we have expressed the solution subsequently according to 

p = Pix;9), (A.3) 

where P{x; 9) is an appropriate function of its arguments. At each instant of time 9 equation 
]3)) defines two curves p(+)(a;; 9) and P(-)(a;; 9) in the phase space such that 



Therefore, if initially the distribution function is given by equation ()A.2|) . then the solution 
of the Vlasov equation (j2.12|) at every subsequent instant of time 9 is represented by the 
expression (j2.14p . 

It is interesting to note that the chain of equations (excluding the last equation) for the 
characteristics ()A.1|) of the Vlasov equation formally coincide with those for the equation 

dP . ^dP . dV ^ f-dU 
o9 ox ox ^ ox 



As a matter of fact, equation ()A.5|) is equivalent to the two equations 

%1 + xpj-f^ = -XX - - Xfp^, (A,6) 

o9 ox ox ^ Ox 

with the initial conditions ()A.4|) . provided U is defined as a solution to the equation ()2.18p . 
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